Regularization of moving boundaries in a Laplacian field 
by a mixed Dirichlet-Neumann boundary condition — exact results 
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The dynamics of ionization fronts that generate a conducting body, are in simplest approximation 
equivalent to viscous fingering without regularization. Going beyond this approximation, we suggest 
that ionization fronts can be modeled by a mixed Dirichlet-Neumann boundary condition. We 
derive exact uniformly propagating solutions of this problem in 2D and construct a single partial 
differential equation governing small perturbations of these solutions. For some parameter value, 
this equation can be solved analytically which shows that the uniformly propagating solution is 
linearly convectively stable. 



Boundaries between two phases that move according 
to the gradient of a Laplacian or diffusive field, occur 
in many fields of the natural sciences and have a long 
and intricate research history Q]; well known examples 
include viscous fingering in Hele-Shaw flow 0, 0] , solid- 
ification fronts in under-cooled melts y|, migration of 
steps Q or electromigration of voids 0, p on the surface 
of layered solids or boundaries of bacterial colonies in an 
external nutrition field "T]. Viscous fingering here takes 
a paradigmatic role as the oldest and most studied prob- 
lem — determining the long time dynamics up to today 
leads to mathematical surprises d, l^, 0, 0| ■ 

A similar moving boundary problem arises in so-called 
streamer discharges 12, 13] that precede sparks and 
lightning. Streamer ionization fronts can be understood 
as moving boundaries separating an ionized phase from a 
non-conducting phase 13, 14, ,15]. The inner front struc- 
ture can be approximated by a boundary condition of 
mixed Dirichlet-Neumann- type, as we will sketch below. 
A similar boundary condition appears in step motion on 
the surface of layered solids when the Schwoebel barrier 
is taken into account 0. Our boundary condition has a 
similar physical effect as the curvature correction in vis- 
cous fingering. We do show here that it indeed stabilizes 
certain uniformly translating shapes. In our analysis be- 
low, we encountered a number of surprises: (z) planar 
fronts are linearly unstable to transversal perturbations 
of arbitrary wave vector < A: < 00, still we find that 
sufficiently curved fronts are linearly convectively stable; 
(m) a simple explicit uniformly translating solution can 
always be found; {Hi) linear perturbations of these so- 
lutions can be reformulated in terms of a single partial 
differential equation, iiv) if the solution of this pertur- 
bation problem is Taylor expanded and truncated at any 
finite order, the eigensolutions seem to be purely oscil- 
lating, {v) however, for a particular parameter value, the 
linear perturbation theory has an explicit analytical so- 
lution that shows that there are no oscillations. Rather, 
perturbations might grow for some time, while they si- 
multaneously are convected to the back where they dis- 



appear. Only a shift of the original shape remains for 
t — > 00. To our knowledge, this is the first explicit so- 
lution showing the convective stabilization of a curved 
front according to the concept of Zeldovich 0] . 

In fact, the interfacial dynamics with our boundary 
condition can be addressed by explicit analysis much 
further than the classical viscous fingering problem. It 
therefore might contribute not only to the understand- 
ing of ionization fronts, but also shed new light on other 
moving boundary problems like the classical viscous fin- 
gering problem. 

A simple moving boundary approximation for a 
streamer ionization front was suggested by Lozansky and 
Firsov 1"^: The front penetrates into a non- ionized and 
electrically neutral region (indicated with a "*") with a ve- 
locity determined by the local electric field E"*" = — Vip"*": 







(1) 



where n is the local normal on the boundary. Approxi- 
mating the interior ionized region as ideally conducting 



ip = const.. 



(2) 



and the electric potential as continuous across the ion- 
ization boundary y)"*" = ip~ , one arrives at the Lozansky- 
Firsov interfacial model. This model was suggested in 
|l3| to explain streamer branching, and it was explicitly 
analyzed in |l,'j ]. Replacing the electric potential by 
the pressure field p, one finds the non-regularized motion 
of viscous fingers in a Hele-Shaw cell. The model gener- 
ically leads within finite time to the formation of cusps, 
i.e., of locations on the interface with vanishing radius of 
curvature [l8| . 

We here propose to replace the boundary condition 



Oby 



Lp — ip = e n • 



(3) 



to suppress these unphysical cusps. Here the length scale 
e characterizes the width of the ionization front where 
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the ionization increases and the electric field decreases. 
It determines the jump ip'^ — of the electric potential 
across the boundary for given field V(/3+ ahead of the 
front. 

The classical boundary condition for viscous fingering 
is — iy9^ = 7K where k, is the local curvature of the 
moving interface, and 7 is surface tension. In contrast, 
the boundary condition Q does not involve front curva- 
ture, but can be derived from planar ionization fronts, 
more precisely from a minimal set of partial differential 
equations for electron and ion densities and their cou- 
pling to the electric field • The formal derivation will 
be given elsewhere. Here we note that ignoring electron 
diffusion [D = 0) as in the planar uniformly trans- 
lating front solutions of the p.d.e.s always yield a relation 
93+ — (p^ — F(r\- V(/5+). For large field i?+ = n • Viy9+ 
ahead of the front, the function F becomes linear, and 
the boundary condition Q results. 

This boundary condition has a similar physical effect as 
the curvature correction in viscous fingering: high local 
fields ahead of the front decrease due to the change of (p+ 
on the boundary, and the interface moves slower than 
an equipotential interface (where = const). While 
the boundary condition of viscous fingering suppresses 
high interfacial curvatures that can lead to high fields, 
the boundary condition ^ suppresses high fields that 
frequently are due to high local curvatures. This physical 
consideration has motivated our present study whether 
the boundary condition Q also regularizes the interfacial 
motion. 

The minimal p.d.e. model for streamer fronts with D = 
leads to a dispersion relation with asymptotes 



be written as 



vk for fc <^ 1/e 
v/e for k :s> 1/e 



(4) 



for linear transversal perturbations 6**^^+** of planar in- 
terfaces . It is important to check whether the moving 
boundary approximation (0-© reproduces this behav- 
ior. Indeed, analyzing planar interfaces we find s{k) — 
vk/ (1 + ek) in full agreement with as we will show in 
detail elsewhere. This strongly suggests that the inter- 
facial model captures the correct physics. It shows fur- 
thermore, that planar fronts are linearly unstable against 
any wave- vector k for all e. 

We now restrict the analysis to the two-dimensional 
version of the model and to arbitrary closed streamer 
shapes in an electric field that becomes homogeneous 

ip{x, y) — » —Eqx far from the ionized body. (5) 

The problem is treated with conformal mapping methods 
[lij : The exterior of the streamer where V^^s^ = can 
be mapped onto the interior of the unit circle. Parame- 
terizing the original space with z = x+iy and the interior 
of the unit disk with uj, the position of the streamer can 



^ oc 

z = X + iy ^ Muj) ^ ^ y2 ak{t) u;'' , (6) 



fc=-i 



where ht{uj) is analytical on the unit disk with a single 
zero at w = and therefore has the Laurent expansion 
given on the right. The boundary of the ionized body 



a e [0,27r[, 



(7) 



is parametrized by the angle a. 

The potential (p^ is a harmonic function due to l|T]l. 
therefore one can find a complex potential ^{z) = ip^+itp 
that is analytic. Its asymptote is ^{z) —Eqz for \z\ 
00 according to (0. For the complex potential ^{uj), this 
means that 



^u) = $(/t(c^)) 



/ 1 °° 
-Soa-i(t) - + V 

\ LO — ' 



Cfc(t) 



(8) 

where the pole oc l/w stems from the constant far field 
Eq, and the remainder is a Taylor expansion that ac- 
counts for the analyticity of $. The boundary motion 
f „ = n ■ Vip ^ is rewritten as 

Re [id^f: dtft] = Re [-z9„$(e^")] . (9) 

The boundary condition Q takes the form 



Re 



= e Re 



(10) 



The moving boundary problem is now reformulated as 
Eqs. and (^Ujl together with the ansatze © and © 
for ft and $. 

For the unregularized problem (where e = 0), it is well 
known that all ellipses with a main axis oriented parallel 
to the external field are uniformly translating solutions: 
for principal radii r^^y = a_i ± oi, they propagate with 
velocity v = —Enlr ^, + rv)/ r,i, while the potential is $ = 

Eoa^i{t)iu;-l/iu)mm- 

For a moving boundary problem with regularization, 

there are only rare cases where analytical solutions can 

be given, and frequently they are given only implicitly 

[l8lll9ll2o|. For the present problem, however, an explicit 

solution is found for all e > 0: 



-Eoa^i 



vt, dta^i 

1 - e/a_i 
l + e/a-i 



0, 



UJ 



(11) 
(12) 



This solution simply describes a circle z ~ x + iy = 
a-ie^*" + vt with radius a_i that according to prop- 
agates with velocity v — — 2i?o/(l + e/a-i). Note that e 
changes the velocity, but not the shape of the solution. 
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Note further that the multiphcity of uniformly translat- 
ing solutions reduces through regularization in a similar 
way as in viscous fingering, namely from a family of el- 
lipse solutions characterized by two continuous parame- 
ters a_i and oi to a family of circle solutions character- 
ized by only one radius a_i or interface width e. 

The physical problem has two length scales, the inter- 
face width e and the circle radius a_i. In the sequel, we 
set a_i = 1, measuring all lengths relative to the radius 
of the circle. 

Now the question arises whether a uniformly translat- 
ing circle is stable against small perturbations, in partic- 
ular, in view of the linear instability of the planar front 
0). The basic equations l|H)l- inT)|l show a quite compli- 
cated structure, and it is a remarkable feature that linear 
stability analysis of the translating circle H11|I " H12() can 
be reduced to solving a single partial differential equa- 
tion. We write 

fticu) = i+r + Piu,T), T^vt, v=-^il3) 
uj 1 H~- e 



-Eq 



1 1 - e 



u 1 + e 



(14) 



where /3 and (j) are analytical in oj and assumed to be 
small. Eqs. © and H1U|I are expanded to first order in /3 
and (j) about the uniformly translating circle and read 

Re [w drP - ujd^l3] = Re [-ujd^(t)] for lj = e''", (15) 



- \ LO 
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- Re [uj'^d^p] = Re [e ujd^cf) + (}>] . (16) 



By construction, F[ijj) = drP — 9^/3 + du)<j) is analytical 
for \uj\ < 1, and Eq. l(TK|l shows that Re[w F{uj)] = for 
cj| = 1. Furthermore, it is clear that uj F{uj) vanishes 
for cj = 0. Therefore, 



^ W F{UJ) = UJ [drP - d^P + 



(17) 



is valid on the whole closed unit disk. The corresponding 
analysis of Eq. (|16|l yields 



^ ( ti> + — 1 uj'^d^P = eujd^4) + + const. 



(18) 



To this equation the operator ujd^ is applied, and Eq. 
(|17|l is used to eliminate terms containing ujduj4'- As a 
result, we find an equation only for the function /3(w, r): 



C, = -e {1 - uj^) UJ dl - {2 + e - 3e uj^) 



(19) 



+2eujdu; dr + 2{l + e) dr. (20) 

Eq. 1)191) has to be solved for arbitrary initial conditions 
/3{uj, 0) that are analytical in some neighborhood of the 
unit disk. The operator Ce conserves analyticity in time, 
e is a singular perturbation that multiplies the highest 
derivatives d^, and d^dr- 



The case e = is almost trivial, since reduces to 

Co ^2{dr- d^). (21) 
Thus all solutions can be written as 



(i{uO,T)=P{LO + T), 



(22) 



where /3(C) is any function analytic in a neighborhood of 
the unit disk |C| < 1. The time evolution just amounts to 
a translation along the strip — 1 < Re C < cx), |Im CI < 1- 
Any singularity of /3 at some finite point C on the strip will 
lead to a breakdown of perturbation theory within finite 
time; this is the generic behavior as found previously in 
the full nonlinear analysis of this unregularized problem. 
Of course, there also exist solutions that stay bounded 
for all times. 

A different perspective on e = is that the Richardson 
moments are an infinite sequence of conserved quantities 
|l8| . A reflection of this property is that any polynomial 
(3{uj,t) = '^^^Qbk{T)uj'' for any N with an appropriate 
choice of the time dependent functions 6^ (r) is an exact 
solution for all times r > 0, for linear perturbation theory 
(|19|l as well as for the full nonlinear problem J^j, i.e., 
any truncation of the Laurent series JBJ leads to exact 
solutions. 

This observation suggests that an expansion in powers 
of a; is a natural ansatz also for nonvanishing (but small) 
e. Taking as initial condition some polynomial of order 
N, one finds from the form (|20|l of that higher modes 
uj'', k > N are generated dynamically — similarly to the 
daughter singularities in regularized viscous fingers 
When the expansion in uj is truncated at some arbitrary 
N' , it can be shown that the problem for any truncation 
N' and for any value e > has purely imaginary tempo- 
ral eigenvalues. One would therefore expect all eigensolu- 
tions for e > to be purely oscillating in time. However, 
this behavior seems inconsistent with our exact solution 
for e = 1. 

For e = \ it turns out that the operator factorizes 



2dr - (1 - c^')a. 



2 + ujd^ 



(23) 



which allows us to construct the general solution. We 
introduce the function 

g(uj,T) = [2^ujd^]P{uj,T), (24) 

that obeys the equation 

[2dr~{l~uj'')d^]g{uj,T) = Q. (25) 

The general solution of this equation reads 

ff(^,T)=G(^^^^), T-tanh^. (26) 

The function G is derived from the initial condition as 
G{uj)=g{uj,Q)^[2 + ujd^](3{uj,Q)- (27) 
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find the asymptotic beiiavior 




G(l) 



(30) 



FIG. 1: Temporal evolution of a perturbed circle ft{(jj) — r = 
1/uj + I3{ijj,t) moving in the positive a;-direction, according 
to Eqs. H12^ and II26|I - I|28|I . The constant motion in time 
T is subtracted. The initial perturbation is a Fourier mode 
Piuj, 0) = -0.5 uj''/{k + 2) with k = 10. The evolution during 
times < r < oo corresponding to < T < 1 is shown for 
time steps T = 0, 0.05, 0.1, . . ., 0.85, 0.9, 0.95, 0.97, 0.99, 1.0. 



independent of to for any initial condition. Therefore 
asymptotically, the perturbation just shifts the basic cir- 
cular solution without change of shape. Indeed, it is eas- 
ily checked that any pronounced structure of the initial 
perturbation that is not located right at the top at w = 1, 
is convected with increasing time toward lo = —1 where 
it vanishes. This is an outflow of the simple dynamics of 
G{Q as pointed out above. Fig. 1 illustrates this behav- 
ior. 

To summarize, we have found that the boundary condi- 
tion ||2Jl at least for e = 1 regularizes our problem in the 
sense that an infinitesimal perturbation of a uniformly 
translating circle stays infinitesimal for all times and van- 
ishes asymptotically for r — cxd up to an infinitesimal 
shift of the complete circle. This statement is based on 
an exact analytical solution for an arbitrary initial per- 
turbation. At the present stage, we have indications that 
this behavior of infinitesimal perturbations is generic for 
e > 0, while the solution is unstable for e = 0. Fur- 
thermore, we expect that the convection of perturbations 
to the back of the structure applies similarly for other 
shapes like fingers. When applying the present calcula- 
tion to streamers, we in fact have to assume this to be 
true, since streamers are typically not closed bodies, but 
rather the tips of ionized channels. Finally, the behavior 
of finite perturbations and their nonlinear analysis will 
require future investigations. 



hence it is analytical in a neighborhood of the unit disk. 
Finally, Eq. (j^^l is solved by 



/5(c^,t) 



dx 



G 



T 



Now the one parameter family of mappings 

LJ + T 



w — > Ct(w) = 



- 1< T < 1, 



(28) 



(29) 



forms a subgroup of the automorphisms of the unit disk. 
Thus on the level of G'(C), the dynamics amounts to a 
conformal mapping of the unit disk < 1 onto itself. 
This dynamics is somewhat distorted by the additional 
integration leading to /3(w,r), but it is easily seen 
that P{uj, t) and d^l3{uj, r) are bounded uniformly in r for 
\uj\ < 1. Hence, contrary to the unregularized problem 
for e = 0, only perturbations contribute that are bounded 
for all times. Hence an infinitesimal perturbation can 
never form cusps. Furthermore, the mapping lo Ct('-^) 
has fixed points uj — ±1; and for r — > oo, i.e., T — > 1, 
it degenerates to Ci('^) = 1> provided uj ^ —1. We thus 
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